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^^ . Abstract. Motivated by the fact that cosmological perturbations of inflationary 

r) I quantum origin were born Gaussian, the search for non-Gaussianities in the cosmic 

(-H ' microwave background (CMB) anisotropies is considered as the privileged probe 

Oj. of non-linear physics in the early universe. Cosmic strings are active sources of 



o 



gravitational perturbations and incessantly produce non-Gaussian distortions in the 



%^ • CMB. Even if, on the currently observed angular scales, they can only contribute a 

C/3 , small fraction of the CMB angular power spectrum, cosmic strings could actually be 

^ ' the main source of its non-Gaussianities. In this article, after having reviewed the basic 

cosmological properties of a string network, we present the signatures Nambu-Goto 

%jj I cosmic strings would induce in various observables ranging from the one-point function 

*vi ' of the temperature anisotropies to the bispectrum and trispectrum. It is shown that 

^^ , string imprints are significantly different than those expected from the primordial type 

Qy ' of non-Gaussianity and could therefore be easily distinguished 

in 

o 

O ! PACS numbers: 98.80. Cq, 98.70.Vc 

, , . 1. Motivations 

c^ ■ The origin of cosmic strings dates back to the discovery that cosmological phase 

transitions triggered by the spontaneous breakdown of the fundamental interaction 
symmetries may form topological defects [H [21 13]. Cosmic strings belong to the class 
of line-like topological defects, as opposed to point-like monopoles and the membrane 
shaped domain walls. As shown by Kibble, the appearance of defects in any field 
theory is related to the topology of the vacuum manifold [3]. If the ground state of 
a field theory experiences a spontaneous breakdown from a symmetry group G to a 
subgroup H, Kibble showed that cosmic strings will be formed if the first homotopy 
group iTi{G/H) 7^ J is non-trivial. In other words, if non-contractile loops can be found 
in the manifold Ai = G/H of equivalent vacua. Similarly, the other homotopy groups 
ttq and 7r2 determine the formation of domain walls and monopoles respectively. Once 
formed and cooled, these defects cannot be unfolded, precisely due to their non-trivial 
topological configuration over the vacuum manifold of the theory. This simple statement 
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suggests that cosmic strings, and topological defects in general, are a natural outcome 
of the unification of the fundamental interactions in the context of Cosmology. As 
remnants of unified forces, their discovery would be an incredible opportunity to probe 
extremely high energy physics with "a telescope" . 

In the last thirty years, many works have been devoted to the cosmological 
consequences, signatures and searches for topological defects [H El El [7]. They 
have pushed cosmic strings to the privileged place to be generically compatible with 
observations. Indeed, domain walls and monopoles are prone to suffer from the 
cosmological catastrophe problem: their formation is sufficiently efficient (or their 
annihilation sufficiently inefficient) to either overdose the universe or spoil the Big- 
Bang Nucleosynthesis (BBN) predictions |8l |9]. For domain walls, this implies that 
either they should be extremely light, i.e. formed at an energy scale less than a few MeV, 
or no discrete symmetry should have been broken during the cooling of the universe. 
There is not so much choice for the monopoles: if interactions were unified, monopoles 
would have been formed. The homotopy group of Ti2{G/H) with H containing the 
U{1) of electroweak interactions is indeed non-triviajj. Cosmic inflation was originally 
designed to solve the monopole problem. If a phase of accelerated expansion of the 
universe occurs, then any defects will be diluted enough to no longer have any (dramatic) 
consequences on cosmology [121 1131 EH IIS]- Meanwhile, Inflationary Cosmology solves 
the flatness and homogeneity problem of the standard Big-Bang model, explains the 
origin and spectrum of the Cosmic Microwave Background (CMB) anisotropics, as the 
formation of the large scale structures [161 E]- Inflation provides a priori an easy 
solution to the topological defects problem by diluting them to at most one per Hubble 
radius. However, one has to keep in mind that this mechanism works only if the defects 
were formed before inflation, and even in that case some may survive [18]. This has to 
be the case for monopoles and heavy walls, but not for local strings. On the contrary, 
exhaustive analysis of particle physics motivated inflationary models, embedding the 
Standard Model SU{?>) x SU{2) x f/(l), has shown that strings are generically produced 
at the end of inflation [19j. In this picture, our universe should contain cosmic strings 
whose properties are closely related to those of the inflaton [20l [211 122] . String Theory 
provides an alternative framework to Field Theories: brane inflationary models propose 
that the accelerated expansion of the universe is induced by the motion of branes 
in warped and compact extra-dimensions [231 1211 [251 126]. Inflation ends when two 
branes collide and such a mechanism again triggers the formation of one-dimensional 
cosmological extended objects, dubbed cosmic superstrings [271 [2H1 |29l [30]. These 
objects may be cosmologically stretched fundamental strings or one-dimensional D- 
brane [311 [32]. Although cosmic superstrings are of a different nature than their 
topological analogue, they produce the same gravitational effects and share similar 

X As often with topological defects, sensitivity to the underlying model is such that one can often find 
a counter-example of any result. Both of these statements, on walls and monopoles, can be evaded in 
some particular models or with some amount of fine-tuning, as for instance if cosmic strings can be 
attached to them and catalyse annihilations [l0| [Tl] . 
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cosmological signatures 

Among the expected signatures, cosmic strings induce temperature anisotropies 
in the CMB with an amphtude typically given by GU , where U is the string energy 
per unit lengthl§| and G the Newton constant [36] • For the Grand Unified Theory 
(GUT) energy scale, one has GU ~ 10~^, which precisely corresponds to the observed 
amplitude of the CMB temperature fluctuations [37]. However, the power spectra do 
not match: topological defects are active sources of gravitational perturbations, i.e. 
they produce perturbations all along the universe history, and cannot produce the 
characteristic coherent patterns of the acoustic peaks [381 Ell HOJIITI |12] . Current CMB 
data analyses including a string contribution suggest that they can only contribute to at 
most 10% of the overall anisotropies on the observed angular scales [l3l HI]. For Abelian 
cosmic strings (see Sec. [2]), numerical simulations in Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetimes show that this corresponds to an upper two-sigma bound 
GU < 7 X 10~^ [15]. Direct detection searches provide less stringent limits but are 
applicable to all cosmic string models: GU < 4 x 10~^ [IHl 1171 HE]. Detecting cosmic 
strings in the CMB certainly requires one to go further than the power spectrum [^ [50] 
(see however Sec. 14. 5p . In fact, strings induce line-like discontinuities in the CMB 
temperature through the so-called Gott-Kaiser-Stebbins effect, which are intrinsically 
of non-Gaussian nature [5T|[52]. In the inflationary picture, cosmological perturbations 
flnd their origin in the quantum fluctuations of the fleld-metric system, and therefore 
were born generically Gaussian. Non-Gaussianities can nevertheless appear from non- 
linear effects during inflation or from couplings to other flelds (see the other articles in 
this issue). These non-Gaussianities are of the primordial type, i.e. they exist before 
the cosmological perturbations reenter the Hubble radius. On the other hand, cosmic 
strings are a source of non-Gaussianity at all times and, as we will see, produce different 
signals from the CMB point of vie^4^ll• 

In this article, we review the non-Gaussian features a cosmological network of 
cosmic strings produce in the CMB anisotropies. In a flrst section, we briefly scan various 
cosmic string models and emphasize their similarities and differences for cosmology. 
Making observable predictions for cosmic strings faces the problem of understanding 
their cosmological evolution. Not only one has to solve the local dynamics in curved 
space, but as extended objects, cosmic strings follow a globally non-local evolution: the 
fate of one string depends on its interactions with the others. The cosmological evolution 
of a network of cosmic strings is a non-trivial problem which can be overcome by means of 
numerical simulations. These simulations permit an estimation of the various statistical 
properties affecting the observational signatures, such as the number of strings per 
Hubble radius, their shapes, velocities or the loop density distribution. Latest results 
in this area, for the Nambu-Goto (NG) type of cosmic strings, are presented in Sec. [31 

§ To avoid any confusion with Greek tensor indices, we will use the Carter's notations U and T for the 
string energy density and tension [35j . 

II Notice that second order perturbations, being non-linear, actively generate non-Gaussianities but at 
a relatively small amplitude [53l [Ml [55l [56] . 
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Once the statistical properties of a cosmological cosmic strings network are known, 
it is possible to extract meaningful observables depending only on the unique model 
parameter U^. In Sec. HI we recap the expected CMB temperature anisotropics induced 
by cosmic strings, derived from various methods. Particular attention is paid to small 
angle CMB maps which preserve all of the projected statistical information. We then 
derive the cosmic string signals expected in various non-Gaussian estimators ranging 
from the one-point function of the CMB temperature fluctuations to the bispectrum 
and trispectrum. We conclude in Sec. [5] and discuss various non-Gaussian aspects which 
still have to be explored. 

2. Cosmic strings of various origins 

Cosmic strings of cosmological interest can be of various kinds depending on the 
microscopic model they stem from. As mentioned in the introduction, they can either 
be non-trivial stable, or metastable, field configurations or more fundamental objects in 
String Theory. From a gravitational point of view, they all are however line-like energy 
density and pressure distributions. In the following, we briefly review the different 
kinds of string having a cosmological interest and we emphasize their similarities and 
differences. 

2.1. Ahelian vortices 

The simplest example of cosmic string illustrating the Kibble mechanism is the Abelian 
Higgs model. The theory is invariant under a local gauge group U{1) and the Higgs 
potential assumes its standard Mexican hat renormalisable form 

vm = \{m'-vi)\ (1) 

where A is the self-coupling constant and rj^ the vacuum expectation value of the Higgs 
field $. In Minkowski space, the Lagrangian reads 

£h = ^ {D^<^)^ (D^$) - \h^,H^^ - V{^), (2) 

where H^i^ is the field strength tensor associated with the vector gauge boson S^ and 

D^ = d^ + igB^. (3) 

At high enough temperature, loop corrections from the thermal bath restore the U{1) 
symmetry and the effective potential has an overall minimum at I'I'I = [H |57] . 
Starting from high enough temperature, one therefore expects the U{1) symmetry to 
be spontaneously broken during the expansion and cooling of the universe. During the 
phase transition, the Higgs field reaches its new vacuum expectation value $ = r/^e*^. 
At each spacetime location, the phase 9{x^) will have a given value, all of them 
being uncorrelated on distances larger than the typical correlation length of the phase 

% If no currents are flowing along the string, Lorentz invariance implies that the string tension T equals 
the energy density U . 
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Figure 1. The Abelian Higgs potential in the complex plane [5ft (<&), 3 (<&)]. The 
non-trivial phase mapping from the internal space to the physical space (right) leads 
to the formation of a cosmic string. The old vacuum j^j = becomes trapped inside 
the new one |$| = rj^. 



transition. As pointed by Kibble, this is at most the horizon size d]^ oc t although one 
expects it to be much smaller [3l |58l [591 EOl [61]. As a result, there exists closed paths in 
space along which 6 varies from to 2tt (or a multiple of 2tt). Such phase configurations 
necessarily encompass a point at which |$| = (see Fig. [1]): the old vacuum has been 
trapped into a non-trivial configuration of the new vacuum, and this prevents its decay. 
Such a structure is invariant by translations along the third spatial dimension and is 
string shaped. 

Solitonic solutions of the field equations describing a static straight Abelian string 
can easily be computed under the Nielsen-Olesen ansatz. The transverse profile of the 
Higgs and gauge field are assumed to be |62] 

QiQ) 



$ = r]^H{Q)e 
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where (r, 9) stands for a polar coordinate system aligned along the string. The 
dimensionless radial coordinate has been defined hy Q = ruhr where m-h = vA?7v is 
the mass of the Higgs boson. The integer n is the "winding number" and gives the 
number of times the Higgs winds the potential for one rotation around the string. From 
Eq. ([2]), the dimensionless equations of motion read 
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Figure 2. String forming field profiles in the Abelian Higgs model with unity winding 
number. The Higgs field H vanishes in the vortex core and reaches its vacuum 
expectation value within distances given by the Higgs Compton wavelength l/wh. 
The gauge field Q condenses inside the vortex over distances given by l/m^,. 



where m\^ = grj^ is the mass of the vector gauge boson. In Fig. [21 we have represented 
the string solution to these equations in Minkowski spacetime [63l [M] ■ The boundary 
conditions are such that the Higgs field vanishes at the center of the string to reach 
its vacuum expectation value (vev) asymptotically. This typically happens after a 
length scale given by its Compton wavelength l/rrih- Similarly, the gauge field boundary 
conditions are such that it has vanishing derivative in the core and remains finite far 
from the string. As shown in Fig. O it actually condenses inside the string with a spatial 
extension roughly equals to l/wib- 

The energy content of such a string is given by the stress tensor stemming from the 
Lagrangian of Eq. ([2]). Along the string worldsheet, 



rptt rpz 



H 



(d^Hf + ^ 
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+ - 
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4 



(6) 



r 4 g- Q" 

which are the only two components which do not vanish after an integration over the 
transverse coordinates. Integrating the temporal part gives the string energy per unit 
length f/, whereas the longitudinal component gives — T. One finally gets 

\J = T = C{\lg^)^l (7) 

where C{\/ g"^) is an order unity function at fixed winding number. Increasing the 
winding number centrifuges the energy density around the core such that U is changed 
in a more complex way [S]. This immediately shows that cosmic strings generically 
carry an energy density and tension of the order of the symmetry breaking energy 
scale ^ ~ r]^. Notice that along the string direction the pressure P^ '- 



-T 



-U is 
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negative, and we are in presence of a "cosmological constant wire" , as one may expect 
from a Lorentz invariant vacuum object. Consequently, the trace of the stress tensor 
Vfj.!^'^'^" vanishes and cosmic strings do not induce any Newtonian gravitational potential. 
Together with the so-called cosmological scaling behaviour (see below), this is the mere 
reason why they remain cosmologically acceptable. They do however induce dynamical 
gravitational effects, the metric far from the string core being Minkowski with a missing 
angle |65] (see Sec. H]). 

The Abelian string model is intensively used in the literature to explore the string 
forming phase transition and string interactions [661 EZ]- From the Kibble's argument, 
one expects the phase of the Higgs field to be random and the resulting string path should 
be a self-avoiding random walk with a given correlation length [68]. Performing lattice 
simulations allows to probe in more details the string forming mechanism and gives a 
more accurate picture of a cosmic string network just after its formation [691 EOl UH E2] • 
Abelian Higgs simulations are also used to compute the cosmological evolution of such 
a network [731 [H EH] (see Sec. Ej). 

2.2. Other flux tubes 

Global string. The Abelian string provides an explicit example of the formation of line- 
like topological defects by the spontaneous breakdown of a gauged symmetry. Breaking 
a U{1) global symmetry can also produce topological defects, the so-called global strings. 
However, in the absence of gauge fields, one can show that global cosmic strings exhibit 
long range interactions and Goldstone radiation [76] . Their dynamics can however mimic 
local strings and being cosmologically acceptable in some regime [771 EHl [79] . 

Non- Abelian string. If the broken symmetry group G is non- Abelian, the cosmic strings 
formed during the phase transition exhibit new properties compared to the U{1) kind |5j. 
In particular, the mapping of the Higgs field to the real space can be made along 
different broken generators $a = r]^ exp{iTa6) . This implies that different type of non- 
Abelian strings may be formed and will interact with each others according to their 
respective windings. The classic example being the appearance of a [Ta,Tf,]-string from 
the crossing between a T^-string and T^-string [80]. In the cosmological framework, new 
strings can potentially be formed at each interaction leading to a frustrated intricated 
configuration [SH [82l [831 [H] . Such an outcome depends on the underlying non- Abelian 
gauge group and Abelian string-like evolution can also be recovered, as for instance in 
the U{N) models [851 [Ml [HIl [HHl [HS] • 

Semi-local string. String-shaped energy density distribution can also appear even if the 
vacuum manifold is simply connected. The non-trivial topology argument is indeed only 
a sufficient condition of defect appearance. The electroweak symmetry breaking scheme 
enters this class, although the first homotopy group is trivial, semi-local strings can be 
formed [901 [HI] • A simple description of these strings can be obtained by replacing the 
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Higgs field in the Abelian model by a doublet in a SU{2) global representation 
These flux tube conflgurations are stabilised because they can be energetically favoured 
for some values of the model parameters, typically for m\, > ruh [93]. Let us notice 
that the currently measured electroweak model parameters do not support stable 
vortices [92]. Similar vortices could also be formed during a chiral symmetry breaking 
phase transition [911 IHSl |96l |971 |98l [99] or within strong interaction models |100[ 1101] . 

K- and DBI-string. These are another extensions of the Abelian Higgs string for which 
the scalar and gauge fleld kinetic terms are non-canonical, or of the Dirac-Born-Infeld 
form |102[ I103[ I104[ 1105] . These strings essentially differ from their Abelian counterparts 
when the gradient terms are non-vanishing, i.e. in the core. 

Current carrying string. In minimal extensions of the Abelian Higgs model, one may 
couple extra-scalar flelds to the string forming Higgs fleld. As shown by Witten, 
this can lead to the condensation of the extra-scalar over the string core |106] . The 
resulting string is carrying a current that breaks the longitudinal Lorentz invariance: 
U and T are no longer degenerate and the string dynamics is affected [351 ESI 1107] . 
One of the most important consequence of these currents is the potential appearance 
of centrifugally supported loops. If stable, these so-called vortons could efficiently 
populate the universe and avoiding the overclosure gives strong constraints on the 
cosmic string energy scale |108[ 1109] . A similar mechanism works for the fermionic 
flelds which are Yukawa coupled to the string forming Higgs fleld. They generically 
produce currents along the string with a discrete mass spectrum, in a way similar to 
the photon propagation in waveguides |110] . Unless the massive propagation modes are 
not excited, the resulting loops are however expected to be unstable |lllt 1112] . 

2.3. Cosmic superstrings 

Cosmic superstrings are fundamental line-shaped objects that can be formed at the end 
of brane-inflation (see Refs. [321I331IM1 111311114] for reviews). The idea that fundamental 
quantum strings can be stretched to cosmological distances has been mentioned by 
Witten |115] . If stable, one would expect fundamental strings to be at an energy scale 
close the String Theory scale, i.e. close to the Planck mass, and this is trivially ruled 
out by observations. In addition, current CMB constraints tell us that the energy scale 
of inflation is at most the GUT scale |116] , implying that strings formed at an higher 
energy would have been diluted anyway. The situation changed with the discovery that 
inflation within String Theory could be a geometrical phenomena induced by the motion 
of a brane moving in a warped throat, somewhere in the compact manifold of the extra- 
dimensions [24j. In the KKLMMT model [261 1117] . the inflaton is a scalar degree of 
freedom associated with the position r of a D3-brane in a warped throat. Within a 
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ten-dimensional super-gravity ansatz for the metric, in the type IIB String Theory, 

ds^ = -^=g^^dx^dx'' + ^/h{r) (dr^ + rM^s^) , (8) 

^h{r) 

the throat is described by the warping function h{r) (exphcitly, it can be the Klebanov- 
Strassler conifold |118] ). In this system, accelerated expansion of our universe comes 
from the interaction of this brane with an anti D3-brane sitting at the bottom of the 
throat tq. Current CMB data suggest that inflation preferentially ends by violation 
of the slow-roll conditions whereas the system continues to evolve till the two branes 
collide [119J . The brane interactions at that stage require String Theory calculations and 
are expected to trigger a reheating era accompanied by a copious production of various 
Dl-branes and fundamental F-strings [271 IMl l29l EQ] . Since the brane annihilation takes 
place at the bottom of the throat, due to the warped metric, the cosmic superstring 
tensions measured by an exterior four-dimensional observer are redshifted by a factor 
/i^/^(ro). The resulting effect is to significantly lower the string tension down to 
acceptable values. In fact, the stability of the produced F-strings and D-strings require 
additional constraints on the model parameters and the spectrum of superstring tensions 
depends on the underlying scenario [30]. For instance, in the KKLMMT model, one 
expects 10"^° < GU < 10"^ [T20] . 

Cosmic superstrings differ from the Abelian strings in various aspects. In addition 
to the coexistence of two different types, they can form bound states of p F-strings and 
q D-strings. The tension of these (p, g)-strings depends on p, g, the binding energy 
but also on their configuration in the throat |12H 1122] . In fact, many of (p, g)-string 
properties mimic the non- Abelian type of topological vortex, as the existence of bound 
states and Y-junctions |123l 11241 11251 [126] . Such similarities have actually been used to 
probe the properties of the cosmic superstrings through the more tractable framework 
of field theory [EH |l28l [I29l [I30] . 

2.4- Infinitely thin strings 

These are the one-dimensional version of the relativistic point particles. Following Carter 
macroscopic covariant approach [35l 11071 11311 1132] , string events can be localised in the 
four-dimensional spacetime by the so-called embedding functions x^ = X'^{^°'), where 
^^ and ^^ are a timelike and spacelike internal coordinate of the string worldsheet. 
Denoting by g the four-dimensional metric tensor, one can define the two-dimensional 
induced metric 

dX^" dX" 

such that the infinitesimal interval between two events reduces to ds^ = jabd^'^d^'' ■ From 
its inverse, on can define the first fundamental tensor 



7ab = ^Mi^i^T^^^TT ' (9) 
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which is nothing but a projector over the string worldsheet. Similarly, -L^i,= g^y — q^^ 
is an orthogonal projector and they verify 

q%q'. = q^. l-%l-',=l-\, q%l-'u=^- (11) 

Variations of the first fundamental form are encoded in the second fundamental tensor 

K./^q^V.q"^, (12) 

where a bar quantity stands for the projection of its four-dimensional analogue over the 
worldsheet, i.e. V^ = Q'^^Vq,. Integrability imposes K, -f = and, by construction, the 
second fundamental form is, respectively, tangent and orthogonal to the worldsheet on 
its first and last indices. As a result, contracting the first two tangential indices gives a 
purely orthogonal vector which measures the string extrinsic curvature |132] 

RP = K^^P = V^g'^". (13) 

The energy content of a spacetime two-dimensional surface can be characterised by 
its internal stress energy tensor. Similarly to the cosmological perfect fluid, one may 
consider a string whose stress-energy tensor is diagonal in a preferred basis. Positivity 
of the energy conditions ensures that the timelike eigenvalue U > 0, while the spacelike 
eigenvalue T should verify \T\ < U |133] . In this frame, U represents the energy per 
unit length of the string and T the string tension. Denoting by m^ and v'^ the respective 
timelike and spacelike orthonormal eigenvectors, one has 

f^u ^ jj^^^u _ j.^^^u = {ij - T)u^u'' - Tq^\ (14) 

where 

u'^Ua = -1, v'^Va = 1, u'^Va = 0, qf"" = -u^u" + w^w^ (15) 

In the absence of external forces, reparametrisation invariance of the string worldsheet 
ensures the stress-energy pseudo-conservation from Noether's theorem [134] 

Vpf "" = 0. (16) 

As for a cosmological fluid, these equations are not sufficient to close the equations 
of motion for the string. One has to supplement them by an equation of state of the 
fluid under scrutiny. The simplest case is the so-called barotropic model for which 
the equation of state is the relation U{T). One can then introduce the two Legendre 
conjugated parameters 

/• df/ , /■ dT 

such that U — T = Jiv. Clearly, v plays the role of a number density and its Legendre 
conjugated parameter p, will therefore be a chemical potential, i.e. an effective mass 
carried per unit number density. Defining their respective worldsheet current density 
by 

Ji'' = Jxu^, u^ = uuP, (18) 
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one can rewrite Eq. fll6p as 

Wpf"^ = fi^Wy + uP\/pJi, - u^V^fip - TK, = 0. (19) 

Contracting Eq. ( TT9l) with v" ensures the current conservation along the string 

Vpv" = 0, (20) 

while its projection onto the worldsheet gives the momentum transport law 

g>^V[p^.] = 0. (21) 

Finally, the orthogonal projection of Eq. ( iT9l) reduces to 



,,. .^., £!V^ .X^, (^ - l) ,', (22) 

where the string acceleration ii^ stands for 

ii" = u'^Vau". (23) 

As should be clear from Eq. (TT4|) . the barotropic equation of state breaks Lorentz 
invariance along the string for U ^ T. In fact, it describes a wide class of elastic string 
models |135[ll36[ll37j . and as suggested by Eq. ( 120|) . the scalar current carrying cosmic 
strings |138] . Conversely, imposing Lorentz invariance along the worldsheet reduces the 
equation of state to the trivial form U = T [see Eq. ([H])], which is also the relation 
found for the Abelian Higgs string. This infinitely thin string is the Nambu-Goto (NG) 
string and does not possess any internal structure |139] . The associated equations of 
motion are purely geometrical and do not depend on U. From Eq. f l22|) . they reduce to 
the vanishing of the extrinsic curvature vector, i.e. 

RP = 0, (24) 

which can be rewritten in a coordinate dependant way by using Eq. (fT3|) 



K>' = -^da iV^r'dbX^') + Vl^y'daX^d^X". (25) 

The connections F^ are for the background spacetime of metric g^^ while 7 is the 
determinant of the induced metric. These equations can also be recovered from the 
usual NG action with an explicit coordinate system |5] 

S=-u[d^^./^. (26) 



3. Cosmological evolution of Nambu— Goto strings 

The previous section shows that the equations of motion of an isolated string depend 
on the underlying microscopic model. The type of string is more determinant when 
two strings interact: cosmic superstrings may form bound states, while non- Abelian 
vortices may weave new vortices from each of their interaction points. Understanding 
the cosmological evolution of a string network requires one to solve both the local 
equations of motion for each string and the outcome of their interactions when 
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they meet. Moreover, the evolution of a system of strings starts from an initial 
configuration which should describe the network configuration just after its formation. 
Numerical simulations have been used to overcome some of the above-mentioned 
difficulties and, up to now, Friedmann-Lemaitre-Robertson-Walker (FLRW) network 
simulations have only been performed with Nambu-Goto strings, Abelian Higgs strings 
and semi-local strings ^M, Hill ^M, [ffl [731 M, EHl IHSl [UHl [751 ^M\. up to some 
variations [7H1 |79l 11291 1148] . As a result, extrapolating the following results to other 
types of string should be made with caution. On the bright side, Eq. ( l22l) suggests that 
as long as the string acceleration remains small compared to T/{U — T), one expects the 
equations of motion of the string to be close to the NG case (up to the eventual vortons 
appearance). In the following, we describe the results obtained for NG strings. Some 
differences exist with the results obtained in the Abelian Higgs string simulations. 

Before entering into details, let us summarize two fundamental properties these 
simulations have revealed. The first is that a cosmic string network avoids cosmological 
domination by evacuating most of its excess energy through some complex mechanisms, 
which typically result in transferring energy between the horizon-sized distances and the 
smaller length scales. For NG simulations, this is the formation of cosmic string loops 
whereas in Abelian Higgs simulation boson radiation is involved. The second property is 
that the infiuence of the initial conditions is expected to disappear on the length scales 
of astrophysical interests. A network of cosmic strings relaxes towards a cosmological 
attractor which depends only on the expansion rate: this is the so-called scaling regime. 

3.1. Dynamics 

The equations of motion for NG strings are the vanishing of the extrinsic curvature 
vector K^^ = 0. In a fiat FLRW background, 

ds^ = a\r]) (-dr/2 + dijdx'dx^) , (27) 

Eq. fl25|) can be simplified with the transverse gauge fixing conditions 



u 



dX'' dX , , 

with the notation t = ^ and o = ^ for the timelike and spacelike string coordinates. 
Such a choice of coordinates refiects the property that a NG string is Lorentz invariant 
along the worldsheet: there is no physical longitudinal component of the string velocity. 
In this gauge, the equations of motion read 
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where a "dot" and a "prime" stand respectively for differentiation with respect to r and 
a. We have also defined the quantity 



X^ 



(30) 



The conditions in Eq. ( 128|) do not completely fix the coordinate degrees of freedom and 
one can supplement them with the so-called temporal gauge fixing which identifies the 
timelike coordinate with the background time at the string event: r = X^ = rj. In the 
transverse temporal gauge Eq. ( I28l) reads X ■ X = 0, while Eq. ( 129|) simplifies to 

= 0, e + 21-LeX =0, (31) 




with 



' 2 

X 



Ti 



X 



(32) 



and H, is the conformal Hubble parameter. The vector symbols being understood as 
three-dimensional spatial vectors. Numerically, it is much more convenient to solve 
an equivalent set of equations found by Bennett and Bouchet |142j . Defining the new 
vectors p and q as 

p(r,M) = ^-X, g(r,t;) = ^ + X, (33) 

evaluated at the new coordinates u = J edcr — r and v = J eda + r, the equations of 
motion (13T]) can be recast into 

dp dq 

— = -H[g-p(p-g)], ^ = -H[p-g(p-g)], (34) 

£ = -n{l-p-q). (35) 

As an illustrative example, these equations have an exact solution in Minkowski space. 
Taking "H = 0, one immediately gets e = 1 (up to a normalisation constant), f){u) and 
q{v) are constant over the characteristics u = a — r and v = a + t. Inverting Eq. fl33|) 
gives 

i(r, a) = ^ [f){a + r) + q{a - r)] , (36) 

which describes the propagation of left and right moving string deformations at the 
speed of light. In the FLRW background, these modes are no longer free moving, but 
interact through the Hubble term [see Eq. (1M1) ]. Solving these equations gives the 
X^(r, a) for each strings but does not predict what happens when two strings collides. 
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Figure 3. Intercommuting strings exclianging tlieir partners (left). On tire right, the 
same mechanism produces a loop from a self-intercommutation. The arrows represent 
velocity vectors. 



3.2. Collisions 

In the infinitely thin approach, the outcome of a NG string intersection event cannot be 
predicted. Stress tensor conservation equations require either that the two strings pass 
through each others, or they intercommute as sketched in Fig. |3l The outcome of a string 
colhsion process can only be addressed within the framework of a microscopic model. 
Numerical simulations of interactions have been performed for a variety of models, and 
in particular for the Abelian Higgs string as represented in Fig. |H In this case, unless 
the relative string velocity is close to unity |149[|150] . or the strings are almost parallel, 
string intercommutation generically occurs |15H 1152] . Let us notice that for type I 
Abelian strings (i.e. having m\, > mh), bound states of the two strings can also be 
formed at low velocity |153[ I154[ 1155] . 

The situation is not necessarily the same for the other types of string. As already 
mentioned, non-Abelian strings can weave new strings from their intersection points 
while current carrying cosmic strings intercommute for bosonic carriers |156] . In the case 
of cosmic superstrings, it has been shown that they intercommute with a probability 
depending on the fundamental string coupling, a quantity which can be significantly 
smaller than unity |157] . In the case of (p, g)-string collisions, Y-junctions can be 
formed under some kinematic constraints |124[ I158[ I155[ I159[ I160[ 1161] . Concerning 
NG simulations, string collisions are actually implemented through a phenomenological 
probability Pe of intercommutation at each intersection event. 

3.3. Initial conditions 

Solving the cosmological evolution of a NG string network amounts to solving Eqs. (!34ll 
and fl35|) along each string, finding all of their intersection points and implementing 
an intercommutation, or not, with the probability Pe- The network evolution is now 
uniquely determined once the initial conditions are specified. The simplest way to set 
initial conditions is through the Vachaspati-Vilenkin (VV) algorithm |162] . Motivated 
by the Kibble mechanism, one assumes a f/(l) Higgs field to be uncorrelated above a 
given correlation length i^. A cosmic string will cross a given plane if one can find 
a closed loop along which its phase runs from to a multiple of 2tt. On a discrete 



Cosmic strings and their induced non-Gaussianities 



15 




Figure 4. Numerical simulations of the intersection between two Abelian Higgs 
strings |152| . The Higgs and gauge field profiles are represented as the inner orange 
and outer pink tubes. Generically, Abelian Higgs strings with rrih > rrih intercommutc. 



three-dimensional lattice, of £c-spacing, it is sufficient to approximate U{1) by Z3 and 
randomly choose the phase at each corner from three values 0, 27r/3 and Att/S to decide 
if a string crosses the associated face. Other symmetry breaking schemes and lattice can 
be approximated in a similar way |163[ I164[ I165[ I166[ 1167] . In Fig. [5l we have shown the 
initial string network configuration obtained from the VV algorithm. The string paths 
have been smoothed by replacing the right angles by circles of radius ic- The initial 
network configuration obviously depends on the physical parameter i^, the network 
correlation length. In FLRW spacetime, there is however another physical parameter 
which has to be specified: the distance to the horizon dh- From those, the initial string 
energy density is now uniquely determined. A random transverse velocity field can also 
be added along each string since one does not expect the strings to be initially at rest 
in any realistic setup. At this point, let us mention that the numerical implementation 
of the VV initial conditions introduce two additional purely numerical parameters: the 
size of the periodic box which contains the simulation, usually normalised to unity in 
comoving coordinates, and the discretisation step required to represent a string, usually 
given by A^ppd, the number of points per correlation length. 
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Figure 5. Initial string network configuration from tire Vachaspati-Vilenkin 
algorithm. The box is a small lattice of 25£c for illustration purpose only. 



3.4- Cosmological scaling 

3.4-1 ■ Long strings. By switching on the evolution from the initial network, string 
motion and intersections drastically change the shape of the strings as well as the 
network aspect (see Fig. [6]). Naively, without any collisional process, one would expect 
the string network to dominate the energy density of the universe. In a volume V , 
denoting by £oo the typical correlation length of the network at a given time (initially 
^oo = 4), the number of strings should be roughly given by V/Cl^. The resulting energy 
density should therefore be 



£L V £L 



(37) 



Due to cosmological expansion one has £oo oc a and poo oc 1/a^. As noted by Kibble, 
this domination does not occur due to intercommutation processes which allow the 
formation of loops. In the so-called "one scale model" , Kibble |3] assumes that loops of 
typical size £oo are formed at a rate equals to l^ (for relativistic speeds, one expects one 
intercommutation per string during the time £oo)- As a result, during a time interval 
(5t, the energy density transferred to loops is 



^Pc 



c^rj5tui^. 



(3^ 
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Figure 6. String network configuration in tlie matter era when tlie horizon fills the 
whole simulation box. String loops with a length smaller than the distance to the 
horizon appear in blue. 



From Eq. (157|) , the energy density of strings which are not loops verifies 



dpo 



dt 



-2ifpoo - 



Pc 



(39) 



where £oo is a function of the cosmic time. Defining C{t) = ioo{t)/t, this equation can 
be recast into 

1 dC 1 f2 + 6w 1 

C dt ' 



(40) 



2t V3 + 3w C^ 

where the background cosmological fiuid sourcing the universe expansion has an equation 
of state P = wp. The constant solution C(t) = (3 + 3w)/{2 + 6w) is an attractor for 
which 

U U 

t2 



Poo OC ^2 °^ ^3(l+t«) • 



(41) 



The energy density associated with strings which are not loops "scales" as matter in the 
matter era and radiation in the radiation era. Therefore, it is prevented to dominate over 
the usual cosmological fiuids and cannot overdose the universe. Notice however that 
the total energy density could still dominate the cosmological dynamics if the energy 
density under the form of loops is not evacuated by some extra-mechanism. For NG 
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Figure 7. Energy density of super-horizon sized strings as a funetion of the eonformal 
time (in unit of Ic) in the radiation and matter era. After some relaxation time, 
Poo oc l/d^j the proportionality factor is universal. 



cosmic strings, loops are transformed into radiation due to the emission of gravitational 
waves |168[ I169[ I170[ 1171] . Other types of loops may loose energy by different radiative 
processes, such as particle emission, or even energy leakage into the extra-dimensions 
in the case of cosmic superstrings [1721 1173] . In Fig. |6l we have represented an evolved 
string network at the end of a matter era run. For FLRW simulations within a 
fixed comoving box with periodic boundary conditions, one cannot evolve the system 
indefinitely: at some point, periodic boundaries become causally connected. Usually, 
one stops the run when the distance to the horizon fills the simulation volume, or more 
rigorously half of it. In this figure, one sees that only a few super-horizon strings remain 
(black long strings) whereas the box is also filled with a lot of small loops (blue) and a 
few larger loops having a size typical of distance between two long strings. The latter are 
freshly formed Kibble loops whereas the existence of the small ones cannot be explained 
in the framework of the one scale model |142] . Concerning the energy density associated 
with the super-horizon strings (also called infinite strings), their evolution in the matter 
and radiation era have been plotted in Fig. [71 they "scale" as expected. From Ref . |144] , 
one has 



Pc 



dl 



28.4 ±0.9, 



Pc 



dl 



37.8 ±1.7. 



(42) 



rad 



As the behaviour of the energy density associated with long strings suggests, the time 
evolution drives the string network towards a stable cosmological configuration which 
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Figure 8. Evolution of the energy density distribution of adh-sized loops as a function 
of the conformal time in the radiation and matter era. After an overproduction regime, 
the energy density distribution scales as l/d\, as the long strings. 



does not seem to depend on its initial configuration, at least for the long strings. Fig. [6] 
therefore displays what a cosmological string network should look like inside a horizon 
volume, at any time during the matter era. As is clear from Fig. [71 the relaxation time 
required for the energy density of long strings to reach the attractor is small. Concerning 
the cosmic string loops, their existence and behaviour have been the subject of various 
claims and analytical works [Hll [HSl [TH ^M, ESI [1771 [UHl [179] . In the following, we 
present recent results |144l 1180] showing that the energy density of loops also reaches a 
scaling evolution similar to Eq. f H2|) . 



3.4-2. Loops. As previously mentioned, the small loops observed in NG simulation 
cannot be explained in the framework of the one scale model. These loops find 
their origin from the self-intercommutation of strings on length scales typical of their 
small size. The building of a small scale structure on strings is the outcome of the 
successive intersection events during which new kinks are produced and propagate along 
the intercommuted segments. Correlations between the kinks induce, from Eqs. ( |33|) . 

auto- and cross-correlations between X{ai) and X(o"2) from which small loops can be 
produced |1741 1177] . In Fig. [HI we have plotted the energy density distribution under 
the form of loops with respect to the conformal time during the radiation and matter 
era. The simulation performed is one of the largest up to date: the box contains lOO^c 



Cosmic strings and their induced non-Gaussianities 



20 



'^ 




Matter era 






E 


10- 


Loop scaling regime spreading to small scales : 


lU 
10" 


Scale free distribution Kibble loops = 
X : 


lo' 

10° 
10"' 


still evolving at time r\ 



«™„(^) 10- 



10 «„ 



Figure 9. Loop number density distribution at the end of a (100£c)^ matter era 
numerical simulation. Apart for a few Kibble loops (aoo < oi < 1), all loops of size 
£ = adh with a > amin follow a scaling regime in which their number density is a 
power law. Smaller loops are still in the relaxation regime and will enter their scaling 
regime later (see Fig. [TOl) . The red dashed line is the best power law fit of the scaling 
function S{a) given by Eq. (|44)) . 



whereas the redshift simulation range reaches almost two orders of magnitude. The loop 
energy density distribution dpo/da is defined such that dpo{a) is the energy density 
carried by all loops having a physical length i in the range adh to {a + da)dh. In other 
words, we measure loop size in unit of the horizon lengthcJ. A logarithmic binning in 
a of resolution Aa/a ~ 10~^ has been used in the range [10~^, 10^] to compute these 
quantities. From this plot, it is clear that after an overproduction regime characterized 
by the bump of Fig. [HI the energy density distribution of loops of given size a relaxes 
towards a stationary regime in which it scales as l/rf^- Such an observation implies that, 
once relaxed, the loop number density distribution is of the form 

dn <S{a) 
da adl' 



(43) 



where the "scahng function" S{a) is found to be well fitted by the power laws |144] 
S{a) = Co a~P with 




1 41 +0.08 
'-■^'- -0.07 



0.09 



-0.03 
+0.03 



and 




1 f,r) +0-21 

-■-•"U.o.is 



0.21 



-0.12 
+0.13 



(44) 



rad 



for the matter and radiation era, respectively. The loop number density distribution, 
for the matter era run, has been plotted in Fig. O As Fig. [H] already shows, the loop 
distribution takes more time to reach the scaling regime for the small loops. The 
relaxation bump is all the more so high and long than a is small. In the loop number 

+ This is the relevant physical length scale of the problem. 
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density distribution, this effect appears as a minimal time decreasing value arr^iniji) such 
that the loop distribution is in scaling at a > amin- The redshift range probes by a 
FLRW string simulation is typically /\z ~ 10^, while for strings formed at the GUT 
energy scales, one expects a Az ^ 10^^ at nucleosynthesis. It is clear that, in the 
cosmo logical context, the string network has quite a time to relax: on all of the relevant 
observable length scales the loop distribution should be in scaling, i.e. a^^^ -C 1. Since 
a power law distribution is scale free, one concludes that a cosmologically stable string 
network does not exhibit loops of a particular size: this is not surprising since the 
only length scale involved is the distance to the horizon. These numerical results can 
be analytically recovered in the framework of the Polchinski-Rocha model [177] . The 
expected loop number density distribution have been explicitly derived by Rocha in 
Ref. |180j with a predicted power p = 1.5 for the matter era and p = 1.8 in the radiation 
era. If not due to statistical errors, these small differences may be explained by the 
existence of additional fractal micro-structure along the strings not considered in the 
analytical approach (see also Sec. H]). Of course, the previous statements hold provided 
the other physical effects which are not included in the simulation do not enter the 
game, as gravitational radiation and gravitational back-reaction. The typical length 
scales at which they should play a role is typically a multiple of GU, or some positive 
power of it |181l I182J . As shown in Ref. [180] . gravitational radiation indeed cures the 



energy density divergence that one can extrapolate from Fig. [9] when a — )■ 0. Let us 
finally notice that although the long strings are defined by a > 1, there is also a small 
population of Kibble loops. Their typical size being the horizon-sized correlation length 
of the long string network, they can be defined to be those having a^o < a < 1, where 

1 fU^'^' 



-'oo 



aoo = -f[ — ] . (45) 

4 \pc 



3. 5. Relaxation towards scaling: memory of the initial conditions 

Although of less cosmological interest, the relaxation of the loop energy density 
distribution towards its scaling regime shows interesting properties which could explain 
some of the differences observed between Abelian Higgs field simulation and NG 
simulations. In the left panel of Fig. [TOl we have plotted the loop energy density 
distribution for loops smaller than the ones in scaling. At the end of the numerical 
simulation, these length scales are those having a < amin- This plot shows that the 
formation of the smallest loops is a delayed mechanism which suggests that a cascading 
process takes place from the initial string network configuration. For a ~ 10~^, only 
the increasing part of the relaxation bump appears at the end of the simulation whereas 
the decaying towards scaling is still visible for the larger a. On the right panel, we have 
plotted the loop number density distribution at the end of the simulation, i.e. at the 
time f]/£c = 50 in the left panel. The loop number density distribution deviates from 
the asymptotic scaling distribution on two typical length scales. 

Firstly, an overdensity compared to scaling is situated at Oc = ic/dh oc 1/t. In 
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Figure 10. Loop energy density distribution with respect of the conformal time for 
the smallest loops. By the end of the numerical simulation, these loops have not yet 
reached their scaling behaviour (left panel) . The right panel is the corresponding loop 
number density distribution at the end of the simulation (jy/^c = 50). Compared to 
the asymptotic scaling power law (red dashed), the loops in the relaxation regime 
are preferentially produced at fixed physical lengths given by the initial correlations 
present in the VV network. 



other words, there is an overproduction of loops with a typical size equals to the initial 
correlation length of the string network. Although one expects the system to retain 
some memory of the initial conditions during the relaxation, it may appear surprising 
that, in spite of the expansion of the universe, the physical length scale of these loops 
remains the same. A physical interpretation is that Cc <^ c^h which suggests that, at those 
small length scales, the system decouples from the Hubble flow. More quantitatively, 
this effect can be explained in the context of the three scale models |174j . Under some 
assumptions, the string small scale correlations can indeed sustain a constant physical 
length. 

Then, there is the overall peak of the loop number density distribution at a^ = 
ir/d^ oc 1/t. Most of the loops which are not in scaling have this size at the end 
of the simulation. In fact, one can check that these loops start appearing soon after 
the beginning of the string evolution. This length scale is, again, at a constant physical 
length 4 and is associated with a purely numerical effect |144j . As discussed, a numerical 
string is discretised with iVppd points. The Bennett-Bouchet code at the basis of the 
simulations presented here uses an adaptive griding algorithm meaning that loops of any 
physical size can be formed |142 j. The only restriction is that, at a given time, a loop is 
an object of at least three points. Consequently, when the initial string network starts 
its evolution, loops smaller than i,- = ic x 3/A^ppci cannot be formed. The existence of 
a finite numerical resolution therefore adds some unwanted initial correlations of length 
i^. Notice that this is not a cut-off but indeed an extra-correlation. 

As a basic consequence, one should not trust a NG simulation at those length 
scales. However, the fact that the initial string network violently relaxes towards scaling 
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by emitting loops at the smallest available correlation lengths has still some physical 
significance |183j . What happens if we increase A^ppd — ?■ oo? As discussed in Ref. [144] . 
the larger length scales are not affected and only the overall peak is shifted around 
the new 4—^0. At these length scales, it is clear that using a NG string to describe 
a network of topological defects would break down and a reasonable assumption is to 
assume that the network will now relax by losing energy through the relevant physical 
mechanism available at those ^r ~^ distances. In Abelian Higgs simulation, most of the 
network energy is emitted through field radiation, up to the point that almost no loop 
are observed in the simulations [731 [75]. Abelian simulation suffers from low resolution 
compared to NG ones and this has been a subject of debate to decide whether or not this 
could explain the absence of loops [TH [75] . The above NG results clearly support that 
particle and/or gravitational waves emission is an important mechanism which certainly 
dominate the relaxation regime. However, when the scaling regime progressively takes 
place, from large to small length scales, the loop formation mechanism should become 
dominant. One may speculate that it is not clearly observed in Abelian simulations 
due to its delayed appearancqj, as it is actually the case in NG simulations when we 
compare it to the formation of 4-sized loops. 

4. Cosmic microwave background anisotropies and non-Gaussianities 

At this point, numerical simulations of cosmic string evolution give us the means to 
derive observable predictions. As should be clear from the previous section, some 
structures in the numerical simulations are not supposed to be present after a realistic 
expansion factor of 10^^: these are the structures still in the relaxation regime and such 
that a < amin, at any simulation time. In the following, we will denote by "systematic 
errors" , the uncertainties inherent to the presence of such non-scaling structures when 
deriving observable predictions from NG numerical simulations. 

4.1. Unequal time correlators 

The first method used to derive CMB anisotropies has been introduced in Ref. |184] 
and applied in Refs. ^SM UMl [M [Ml [771 HHS [HH] for global topological defects and 
recently in Ref. |190l 1191] for the Abelian strings. Cosmic strings are active sources of 
gravitational perturbations [391 1193] which means that the equations of motion of their 
induced linear perturbation is of the form 

VV = S, (46) 

where "D is a time differential operator, V the perturbation in the energy density, or 
velocity, etc. . . which is directly related to the CMB temperature anisotropies. Here S 

* A hint for this is the compatibihty of the Abehan Higgs string correlators with the Polchinski-Rocha 
model [75 ; this one explaining the NG loop distribution }18Q] . 
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Figure 11. CMB temperature power spectrum induced by a network of Abelian 
cosmic strings and derived from the unequal time correlators method. Figure extracted 
fromRef. [192]. 



denotes the source terms, i.e. the string stress tensor. From the Green's function Q of 
this equation, one gets, today (at r/o) and in Fourier space 



■no 



V{k,r]o)(x / Qk(r])S(r],k)<lrj. 



(47) 



The two-point correlator reads 

r rvo 
{V\r^„ k)V{vo, ^))^ j j Gliv')Gk{v){S\v', k)S{v, k))drjdv', (48) 

and its determination requires a full-time knowledge of the source term S{ri, k) for each 
mode. Since it is impossible to carry out a simulation over the whole cosmological 
history, the scaling properties of the cosmic string network can be used to analytically 
extrapolate the source terms over the required ranges. As shown in Refs. |187j . as long 
as the cosmic string network is in a scaling regime, the source terms are the stress tensor 
components and assume the form 

(r^.(A;,r/)r;,(fc,r^'))« ^ 



V 






(49) 



Numerical simulations are actually used to determine the scaling functions /, 



fivpa- 



In 



Fig. [m we have represented the temperature anisotropics derived in Ref. |192] using 
such a method from Abelian Higgs string simulations. The current CMB constraint on 
GU comes from this power spectrum: at two-sigma, GU < 7 x 10^^ [194j. 

4-2. Simulated small angle maps 

The previous constraint typically corresponds to a string contribution which cannot 
exceed 10% at the multipole moment i = 10. On current observable angular scales. 
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Figure 12. Stacking of two numerical simulations to compute the small angles CMB 
anisotropics induced by NG strings |199j . The dashed rectangles represent the redshift 
and angular extension of two numerical simulations used to evaluate Eq. (1521) . The 
first starts at last scattering and stops at z ~ 30 while the second starts at 2; ~ 30 and 
ends now, up to a small overlapping to ensure relaxation towards scaling. 



cosmic strings may only be a sub-dominant fraction of the overall CMB anisotropics. 
However, string induced perturbations being non-Gaussian, as opposed to inflationary 
perturbations of quantum origin, one can go further than deriving the two-point 
function. Notice that, in principle, the unequal time correlator approach could be used 
to extrapolate the three- and higher n-point function by using the scaling properties of 
the string network. Another approach is to produce simulated maps of string induced 
CMB anisotropics. Again, we face the problem of the small redshift range probed by 
the numerical simulations. By putting an observer inside the numerical simulation, such 
maps can only include stringy effects up to a finite redshift, typically z ~ 10^. The CMB 
anisotropics computed in this way are therefore only accurate on large angular scales 
but can produce full sky maps |195[ I196[ 11971 1198] . This limitation can be avoided by 
stacking maps from different redshifts, an approach outlined in Ref. [36] and applied in 
Ref. [T99l[2nn] . 

Simulations with the observer outside of the numerical box are not well suited for a 
full-sky map reconstruction, but are perfectly designed for the small angular scales. The 
reason being that cosmic strings are incessantly sourcing the CMB fluctuations since 
the last scattering surface, and contrary to the perturbations of inflationary origin, 
this part cannot be affected by Silk damping. Therefore, at small angular scales, one 
expects the strings' signature in the CMB temperature fluctuations to be dominated by 
their integrated Sachs- Wolfe (ISW) effect from the last scattering surface [200] . In the 
temporal gauge (r = X^ = rj), the NG stress tensor derived from Eq. (126!) reads 



rpfiu 



u 



da eX'^X'' --X^XnS 



x-X). (50) 



-9 J \ ^ 

In the flat sky approximation, well suited for angles typically smaller than the Hubble 
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angular size at the epoch of interest, Hindmarsh has shown that the ISW temperature 
anisotropies induced by NG strings can be simphfied to |20H 1202] 

k J X n x^ 
where Q{x) = 6T(x)/Tcmb- The wave vector k denotes the transverse component of the 
three-dimensional vector k with respect to the line of sight n, whereas, in the temporal 
gauge, M encodes the string stress tensor distortions of the photon temperature and 
reads 

^ --^ (n-X)-X 

u = X-^ ^^. (52) 

1 + h-X 
As can be seen in Eq. f ISip . only the strings that intercept the photon path x^ can 
imprint their signature in the CMB temperature fluctuations. The previous expression 
is nothing but the Gott-Kaiser-Stebbins effect in the temporal gauge [52l I20H 1203] . As 
a result, the knowledge of u, and therefore of the string trajectories X, is only required 
on our past light cone. In the context of string numerical simulations, the trajectories 
of all strings are computed during all of the numerical simulation time. Therefore, to 
compute u, one only needs to determine which parts of the string network intercept our 
past light cone and at what time. 

In order to be able to generate a significant amount of maps using such a method, 
it is more convenient to stack two medium sized (50£c)^ simulations along the lines 
sketched in Fig. [121 The first one starts at the last scattering surface and ends at 
a redshift fixed by the maximum expansion factor achievable in the numerical box. 
For the simulations we performed, initially dh ~ 0.185 (in unit of the fixed comoving 
simulation volume), which corresponds to 1.7 Gpc and a field of view of 6'fov — 7.2° 
(for the current fiducial cosmological parameters [42j). Such a run ends after a 30-fold 
increase in expansion factor, corresponding to a redshift z ~ 36. We then propagate 
the photons perturbed by the first run into a second numerical simulation of the same 
size but starting at Zi ~ 36. For another 30-fold increase in expansion factor, this run 
ends at 2; ~ 0.3. As can be seen in Fig. [121 the second simulation represents a much 
larger real volume than the first one and therefore subtends a greater angle in the sky. 
As a result, only the sub-part of the second run that matches the angle subtended by 
the first simulation is actually used. As we will see later on, the CMB temperature 
maps are weakly sensitive to the string network at low redshifts, simply because there 
are almost no strings intercepting our past light cone in a recent past, which makes 
this technique perfectly acceptable. In practice, each of these numerical simulations 
is started before the redshifts mentioned, in order to give the cosmic string network 
enough time to relax toward its stable cosmological configuration. As discussed in the 
previous section, one has to make sure that the structures (strings and loops) we are 
interested in have indeed reached their scaling behavior during the numerical runs. We 
switch on the photon propagation inside the runs only after making sure all the large 
structures (infinite strings and loops) are in their scaling regime. This can be checked 
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Figure 13. String-induced CMB temperature fluctuations on a 7.2° field of view [44) . 
Because of their cosmological scaling, most of the long strings intercept our past light 
cone close to the last scattering surface. As can be seen in the right image, the edges in 
the temperature patterns can be identified to strings intercepting our past light cone. 
Note that active regions corresponding to string intersection and loop formation events 
lead to the bright spots in these maps. 



by monitoring the evolution of the energy density distributions, and we have chosen to 
start the photons' propagation when all loops larger than a third of the horizon size are 
in scaling. This cutoff is then dynamically pushed to smaller values to include all the 
loops entering the scaling regime at later times. The cutoff time dependence is simply 
the function amin(?/) and can be deduced from the loop distribution relaxation times 
derived in Sec. 13.4.21 The resulting CMB temperature map is displayed in the left panel 
of Fig. [13] whereas the right panel shows the string paths projected onto our past light 
cone. Again, these maps are only representative at small angular scales. On larger 
angles, they represent only the ISW contribution to the overall string anisotropics: for 
instance, Doppler contributions coming from photon decoupling at last scattering are 
dominant around I = 300 — 400 (see Fig. [TTl) . 

The discussion on systematic effects coming from the presence of loops not yet in 
scaling can be found in Ref. |44]. In fact, they have only a small effect. The physical 
reason being that, due to scaling, the long strings are still the main source of CMB 
anisotropics even at (reasonably) small angles. Indeed, there are always roughly ten 
strings per Hubble volume at each time, which means that a patch of 0.8° is at least 
crossed by the ten long strings being there at last scattering, plus a few others from 
lower redshifts. 
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Figure 14. The probability distribution function of CMB temperature fluctuations 
induced by NG cosmic strings. The orange curves quantify systematic errors coming 
from the string simulations by including non-scaling loops, or by removing all loops. 
Deviations from Gaussianity are clearly apparent in the tails of the distribution, as 
well as from the negative skewness. 



4-3. Skewness and kurtosis 

The most basic statistical test that can be performed from a set of small angle CMB 
maps is to plot the one-point function of the temperature anisotropics. As can be 
seen in Fig. [TH the temperature anisotropics induced by cosmic strings are clearly 
non-Gaussian. From a set of 300 independent CMB maps, one finds the mean sample 
skewness to be negative 



9i 



(6-0)3 



cr^ 



-0.22 ±0.12, 



(53) 



where the brackets stand for the mean over different realisations while the bar denotes 
averaging on each map. The variance itself averages to 

(T^ = /(e-e)2\ ~ (150.7 ± 18) {GUf. (54) 

The quoted errors are statistical and refer to the square root of the variance between 
the different realisations. Similarly, the mean kurtosis averages to 



92 



(6 - 0)^ 



o^ 



0.69 ±0.29. 



(55) 



An analytical approach extending these results to cosmic superstrings can be found in 
Ref . |204] . A simple way to look for strings is to search for large (but rare) temperature 
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fluctuations. Deviations from Gaussianity start to be significant, let us say by a factor 
of two, only in the tails when the probability distribution becomes typically lower than 
10-*^. 

4-4- Real space methods 

Strings induce step-like discontinuities in the CMB anisotropics and various methods 
have been designed to probe the non-Gaussianities associated with them. 

Multifractal analysis |205j has the advantage of being directly applicable to the 
time-ordered data retrieved when a CMB telescope scans the sky. As opposed to the 
fractal dimension of a set which measures how sparse it is, the multifractal spectrum 
of a measure defined over a set gives how many and which fractal dimensions there 
are. In the context of cosmic strings, this method has been applied in Ref. |206j on one 
dimensional scans of maps similar to the one in Fig. [131 the measure being defined by 

MO = [0W- 0(^ + 1)]', (56) 

where the integer i labels a point along the scan. The multifractal properties of this 
measure have been shown to be distinctive enough to detect strings, compared to a 
Gaussian signal, but only when the detector resolution is sufficiently good. One may 
wonder, under multifractality, how a non-Gaussian string pattern could be distinguished 
from other non-Gaussian sources. In fact, Eq. ( l56l) consists in taking the gradient of 
the induced CMB fluctuations along the scan. Step-like discontinuities, passed over a 
gradient fllter, become one-dimensional delta functions, and this is a definite string 
feature that can only be altered by the beam experiment. Denoting by a and /3 
the horizontal and vertical angular coordinates, the gradient magnitude |V0| of the 
temperature anisotropics is defined by 



This definition makes it clear that for a finite temperature step, let's say G(a,/3) = 
OoH(a — ao), H being the Heaviside function, the resulting gradient magnitude is a 
Dirac distribution at the string location. In Fig. [151 we have plotted the gradient 
magnitude of the temperature maps of Fig. [131 ^'S well as its convolved version with 
a Gaussian beam typical of the Planck satellite at 217 GHz. With a finite resolution 
beam, the discontinuities are now smoothed. Real space methods applied to string 
are therefore strongly sensitive to the angular resolution. Let us mention that wavelet 
analysis methods have been also explored in this context [207[ I208J or to produce cleaner 
maps [209] . 

Directional gradients, obtained by variations with respect to either a or /3, have been 
discussed in Ref. [203j in the context of Minkowski functionals. They are again found to 
provide a more distinctive non-Gaussian behaviour than the original temperature map 
for the contour length and genus. 
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Figure 15. Normalized gradient magnitude of the string-induced temperature 
anisotropies shown in Fig. [13] (left). A logarithmic scale has been used to enhance 
the contrast by preventing the bright spots from saturating the color scale. The right 
panel is the gradient magnitude obtained after convolution by a Planck-like Gaussian 
beam of resolution 5'. Notice that the color scale is back to linear, most of the bright 
spots being now smoothed by the beam. 



4-5. Temperature power spectrum 

Moving to Fourier space, the small angle CMB maps also permit a determination of the 
power spectrum at large multipoles. In Fig. [161 we have plotted its mean value over the 
different maps as well as the one-sigma statistical error around the mean. The overall 
power at £ = 1000 is [B] 

m + 1) G 



27r 



u{Guy 



{hi 



^=1000 



which is close to the value obtained in Abelian Higgs field simulation (see Fig. [TT]) . This 
is not so surprising since the long strings in both NG and Abelian Higgs simulation have 
a similar scaling evolution, and as explained above, long strings are the main sources of 
CMB anisotropies even at the small angles. The power law tail in Fig. [16] is the direct 
consequence of the presence of strings at all times since the last scattering surface: one 
finds for 1^1 



e{i + i)Ci(x t 



with 



P 



0.889lo:o90 ) 



(59) 



where only the systematic errors have been reported. Such a power law shows that 
cosmic strings have to become the dominant primary source of CMB anisotropies at the 
small angular scales, the fluctuations of inflationary origin being killed by Silk damping 
at those multipoles. In Fig. [TTJ we have plotted the respective contributions of strings 
and adiabatic anisotropies of inflationary origin. The cosmological parameters have 
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Figure 16. Mean angular power spectrum of the string induced CMB anisotropics at 
small angular scales and its one sigma statistical errors (averaged over 300 maps). At 
"small" multipoles I < 500, Doppler contributions from the last scattering surface are 
expected to be significant and this plot gives only the ISW component [44] . 
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Figure 17. CMB temperature anisotropics from various sources, compared to the 
expected string contribution. Due to the Silk damping of primordial perturbations, 
string induced anisotropics always become dominant for the large multipoles. However, 
an unresolved SZ effect may compromise such a clean signature |44j. 



been set to their fiducial values in the Lambda-Cold-Dark-Matter (LCDM) model and 
the string energy density U is compatible with the current bounds. For the current 
upper limit on GU = 7 x 10~^, CMB anisotropics should become dominated by cosmic 
strings at £ > 3000. An unresolved Sunyaev-Zel'dovich (SZ) component may, however, 
compromise such a signal. Nevertheless, string induced anisotropics do not depend on 
the signal frequency whereas the SZ does, and one may hope to disentangle both 
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4-6. Hindmarsh approximation 

As shown by Hindmarsh in Ref . |201j , the power law behaviour of the NG string power 
spectrum at small angles can be analytically recovered. In Fourier space, the power 
spectrum of the string induced CMB anisotropics 9 is defined by 

' Qk.Qk,) = P{h){27i)H{k^ + k2), (60) 



the expression of Gfe being given by Eq. (IHT]) . Remembering that the string positions 
and velocity vectors have to be evaluated on the past light cone, it is more convenient 
to use the so-called light cone gauge. Instead of identifying the timelike worldsheet 
coordinate r with the background time at the string event, one chooses instead to 
identify r = X+ = X° + X^. In this gauge, Eq. (15T!) simplifies to 

- k'^Qk = itkA j daX^(a)e*'=-^("), (61) 

where we have defined 

e = SyrGf/, (62) 

and where the capital indices are two-dimensional. The time parameter r then labels 
the intersections of a set of null hyperplanes with the worldsheet. For our problem, all 
quantities have to be evaluated a,t t = x^ = t] + z. In a field of view of formal area 
A = (27r)^5(0), one can express the power spectrum as 

P^k) = e'^ fdada' (x^(a)X^(a')e^'-[^M-^('^')]) . (63) 

Adding the assumptions that both X"^ and X^ obey Gaussian statistics, all of the 
correlation functions of G^ can now be written in terms of two-point functions only. 
Using the same notation as in Ref. \20T\ . the non- vanishing two-point functions are 

X^(a)X^(a'))=V^l^(a-a'), (64) 

X^(a)X^(a')) = V^M(a-a'), (65) 



X^(a)X^(aO) =^5^^T(a-a'), (66) 

as well as the quantities 

r(a - a') = ([X(a) - X(a')]') = J'^da^ j'^da2T{a^ - a^), (67) 

U{a-a')= (^[X{a)-X{a'))]-X{a')^= T daiM(ai - a'). (68) 

The leading terms are given by |201l I210J 

[Oa— >oo l^aa— J-oo 

n(.) - { ^f ^ - ° , (TO, 

a — !■ oo 
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where we have defined 

e = r'(oo), v^ = lx'), t' = lx'), co = e(i:-x\. (71) 



The correlation length ^ is the projected correlation length on the past light cone, P is 
the mean square projected tangent vector, w^ is the mean square projected velocity and 
Co the correlation between projected velocity and curvature. From these assumptions, 
Eq. (163|) reduces to 



P(fc) = ^-^ /dada' 



V{a - a') + -PU^a - a') 



^^k^r(a~a')/4_ (72) 



When kC, gets large, the terms involving the mixed correlator M can be shown to be 
sub-dominant and only the first term remains: 

P{k) = ^J da^da^V{a^)e-''^^-y\ (73) 

where a± = a ± a' . Denoting by L the total transverse light-cone gauge length of string 
in the box of area A, one gets 

k^P{k)c.e^^^^^-^ (74) 

A t (^kO 

At small angles, the wave number fc^ ~ £{£+ 1) and Eq. (|74l) predicts that £{£ + l)Ci oc 
£~'^. The small difference with Eq. 059 p is suggestive of a cloud of zero- dimensional 
objects along the string worldsheet which may be the signature of small loop production 
in the NG numerical simulations. Let us stress that Eq. (1711) is not "primordial" but 
directly approximate the observed angular power spectrum of the CMB temperature 
anisotropies. 

4-7. Bispectrum 

4-7.1. Analytical approach. The success of Hindmarsh approximation to describe the 
small angular CMB anisotropies power spectrum suggests it can be applied to higher 
n-point functions. In Ref. |210j . this method was used to derive the bispectrum defined 
from the three points function by 

{Qk,Qk,Qks) = Bik,, fc2, k,){2ny6ik, + fc2 + fcg). (75) 

Plugging Eq. fl6T|) into the previous expression gives 



B{k^, k,, fcs) = ^^'^^^^10^ I da,da,das (x^X^X^e^'^''-''^) , (76) 

with X^ = X^{aa), a,b E {1,2, 3}, and ^1+^2+^3 = 0. With the Gaussian assumption, 
the ensemble average of the string observables is lengthy but straightforward and the 
final result reads |210j 

B{ki,k2,k3) 



° t4 ^ ^"2 kjkikl 



kfK23 + k^Ksi + k^Ku 
{fi23f^31 + /tl2K3i + K12K23) 



(77) 
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The quantities i^ab are shorthand for the scalar products Kah = —ka-ki,. In the same way 
as for the power spectrum, this expression directly gives the bispectrum of the CMB 
temperature anisotropics. Its overall dependence varies as 1/k^. Its sign depends on the 
sign of Co defined in Eq. (17T]) . and contrary to what one could naively expect Cq 7^ 0: the 
projected string velocity and curvature vectors are correlated. This can be shown by 
starting again from the equations of motion fl29l) . but this time, in the light cone gauge. 
The equation of motion for X^ gives 

- + 2H(X° + X^) = 0, (78) 

whereas the equation for the transverse components is 

' e'^ V J e da \ 6 da I e 



x + 2n^ix)x---i---\-2n-,ix.x)x = o. (79) 



. 2 

In a FLRW background, assuming that {X ) is constant, and neglecting higher-order 
correlations between "H, X and X, we find 

30) 




ds"^ 

where we have defined ds = eda, and where 1-L is the averaged conformal Hubble 
parameter. Still assuming that the ensemble is approximately Gaussian in X and X /e, 
the right hand side reduces to 

'^"'' .x\=h((x'\(£\-(x.x''^ 



ds'' 

The last term vanishes and the cross correlator Tiv'^P is positive: from Eq. fl7T]) . we 
deduce that cq > 0. It is interesting to notice that cq would vanish in Minkowski 
spacetime, which can be viewed as a consequence of time reversal invariance. The 
existence of a cosmic string bispectrum is the consequence of the breaking of the time 
reversal invariance in a FLRW background. 

An illustrative example is to apply Eq. f l77|) to the isosceles triangle configurations 
in Fourier space such that 

n 

\ki\ = \k2\ = k, Ifcsl = 2A;sin-, (82) 

where 9 denotes the angle between the wavevectors ki and ^2. The isosceles bispectrum 
reads 



v'^L^ 1 l + 4cos^sin2(^/2) 
'¥^'^ sin^^ 



Bueik, 6) = -e^vrco^^^— .3. • (§3) 



Notice that for 9 = 7r/3, we obtain the peculiar case of an equilateral triangle. 
In Fig. [IHI, we have plotted the angle dependency of the isosceles bispectrum. These 
configurations are amplified as 1/9^ in the two flat triangle limits for which either — )■ 
or 6' — )■ vr. Both of these configurations are therefore better suited than the equilateral 
one to characterize the strings. As suggested by the real space searches, the strings can 
produce a strong bispectrum signal only if the detector resolution is sufficiently good. 
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Figure 18. Angular dependency of the isosceles bispectrum as a function of the angle 
9 in between the wavevectors fci and ^2- The particular values 6'oq — 7r/3 corresponds 
to the equilateral configuration and 9o makes the bispectrum vanishing. Notice the 
amplification for fiat triangle configurations at 6* ^- (squeezed) and 6 —>■ tt (collapsed). 



Assuming a beam resolution of 5' means that the 7.2° field of view would contain at 
maximum roughly 80^ Fourier modes. Consequently, the smallest values of 6 achievable 
would be around 6 > 0.03 radians, with only a few modes saturating this bound. 



4- 7.2. Numerical results. The previous analytical results can be compared to the CMB 
temperature bispectrum derived from the simulated maps of Sec. 14.21 Numerically, one 
can use the scale convolution method introduced in Ref. |21H 1212] and applied to the 
string bispectrum in Ref. |210] . This method relies on the choice of unity window 
functions in Fourier space Wu{l) peaked around a particular wavenumber u. Defining 

r AI - 

0«N = y (^0^W;(/)e-'-, (84) 

one can construct an estimator of the three point function in Fourier space by remarking 
that 



dpdqdk 



(2vr)^ 



&p&,&kWk,{p) WkM Wk,{k) {27ryS{p + q + k). 



15) 



For thin enough window functions, Qk remains constant over the window function width 
and we construct our reduced bispectrum estimator as 

1 



Jk-ik2kz 



s, 



k\k2k-i 



&kiix)Qk2ix)Qk3{x)dx ) . 



(86) 



The function S^'^^ is the fiat sky equivalent of the inverse Wigner-3j symbols and reads 



'felfc2fc3 



^,WkMWk2{q)Wk,{\p + q\). 



(87) 
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Figure 19. Left panel: mean value and standard deviation of the squeezed isosceles 
bispectrum [£(^ + l)/(27r)] ' bug for 9 = 0.2 radians. The dashed line is the best 
power law fit. The right panel shows its rescaling by 0^[(.{l. + l)/(27r)]^/^6«6( showing 
the 1/9^ dependency. The spurious plateau (dotted) for the lower multipoles conies 
from a numerical cutoff associated with the window functions and occurs at fmin(^) — 
307r/(e'fov6'), the field of view being 0fov = 7.2°. 



For the window functions such that Wu{k) = 1 for u — w/2 < k < u + w/2, one can 
approximates 

Wuik)c:^wS{k-u), (88) 



for small enough values of w compared to the wavenumber k. In this case, Eq. (jH 
be worked out into 

3 4kik2k3 



^'- - (S 



can 



^9) 



[(A;i + k^y - kl] [kl - (fci - k^r] 

In the left panel of Fig. [111 we have plotted the mean string bispectrum and 
its standard deviation obtained by this method over the 300 string CMB maps. For 
illustration purpose, this plot is for the isosceles configuration having 6 = 0.2 radians. 
The right frame of Fig. [TI?] shows the same mean bispectrum but multiplied by 6^, for 
various small values of 6. As expected from the analytical results, we recover the 1/6^ 
behaviour. The wavenumber dependency also matches with the analytical calculations, 
up to similar slight power differences as we found for the power spectrum. A power law 
fit against the mean numerical estimator gives 

[i{i + l)f^ hue oc £-" with g ~ 2.8 , (90) 



while the overall amplitude can be evaluated around the minimum variance multipole. 
At £ = 5000, one gets 

3 



+ l)/{27,)fHue 



=5000 



;-2.7±1.4) X 10" 



GU 



(91) 



which also matches with Eq. ( 15^ under some crude estimation of the string 
parameters |210] . Finally, as suggested by Fig. [181 the string bispectrum is mostly 
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negative. Integrated over all possible configurations, one recovers the mean negative 
sample skewness of Eq. ( !53|) . thereby explaining its origin as a direct consequence of the 
breaking of the time reversal symmetry in FLRW spacetimes. 

4-8. Trispectrum 

The trispectrum of the string induced CMB temperature anisotropies can be derived in 
a similar way. Starting from the definition of the four-point functions 

[Qk.Qk.Ok^Qk,) = T{k^, k2, ks, fc4)(2vr)2(5(fci + ^2 + ^3 + k^), (92) 

we define the trispectrum a^ 

T{ki,k2,k3,k^) 
_ e ki^k2gk3^,k4^ 

withXf = X^{aa), {a,b) G {1, 2, 3,4} and fci + A;2 + fc3 + fc4 = 0. As shown in Ref. [213] . 
the trispectrum and the higher n-point functions exhibit unfactorisable "fiat directions" 
in the n-dimensional space of the integration variables {c^a}- Physically, it means that 
the leading order part of the (connected) trispectrum is sensitive to the higher orders 
of the correlators in Eqs. ( IMl) to ( 166|) . For the trispectrum, the correlator T{a) has to 
be expanded at next-to-leading order, and following the Polchinski-Rocha model [177i| . 
we assume a non- analytical behaviour for T{a) at small scales 

2X 



d(Tid(T2d(T3da4 (x^XiX^X^e'^ ''-■'''' ) , (93) 



72 ' ^ 



T{a)^t'-c,i^jj . (94) 

In the light-cone gauge, we leave Ci and x ^^ undetermined parameters since they 
cannot be straightforwardly inferred from the numerics performed in the temporal gauge. 
Nevertheless, because the correlation should be smaller as a becomes larger, one should 
have ci > 0. In the temporal gauge, x is directly related to the power law exponent of 
the scaling loop distribution functions through x = 1 — p/2. As we are bound to show, 
the mode dependence of the trispectrum will also be uniquely given by this parameter. 
Once the tangent vector correlator expressed as in Eq. ( l94ll . the integrations in Eq. (|93ll 
can be performed explicitly, except for parallelogram configurations which have to be 
dealt as a special case. After some tedious calculations, an interpolating expression for 
the trispectrum is |213] 

Tik,,k2,k„k,)^e'^^^ (cif)"'^^''^'V(x) 

X 7n (^^' C2Y'A'^A g{k,, fc2, fc3, k,). (95) 

In this equation, f{x) is a number depending only on the parameter x 

fix) = ^r (^-^) [4(2x + l)(x + 1)]^/^'^^'^ , (96) 

jl Notice that our denomination "trispectrum" here contains the unconnected part. This one is however 
non- vanishing only for paraUelogram configurations of the wavevectors. 
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while g{{ka\) is the trispectrum geometrical factor defined by 

g{ki, /C2, fcs, ^4) = ,2,2.2.2 [^ \ ' (97) 

rb-| /Co /Co /C^ 

where 

r2(fci, fc2, fcs, ^4) = -K12 {klkl - 4^Y^^ + C, (98) 

and O stands for cyclic permutations over the indices. The function 7n(a, x) stands for 
the normalised incomplete lower gamma function defined by 



7n(a, x) ^ ^^-Y = T^ / t^-'e-'dt, (99) 

r(a) r(a) Jo 

and, finally, A has been defined by 

A(fci, fc2, fcs, ^4) = TTo X -^^^ -(100) 

{klkj - KJ^) ' + O '*i2/t34 + «;i3«;24 + ^^\^^^2^ 

As an application, the trispectrum over parallelogram configurations is obtained 
when y^ = and the leading term of Eq. f l95|) simplifies to 

To(*..fe.fc,*.)-^^^p|-^. (101) 

where Q now refers to the parallelogram angle. Under the scaling transformation 
ka — !• Afca, the parallelogram trispectrum scales as 

To(Afci, Afc2, Afca, Afc4) = A-%(fci, ^2, fcg, ^4). (102) 

For parallelograms, as already mentioned, the trispectrum also gets a contribution from 
the unconnected part of the four-point function, which is Gaussian and reads 

ro-(fci, fc2, fc3, ^4) = AP{k^)P(k.,)^ O . (103) 

Using Eq. ( 163|) . ones sees that the unconnected part also behaves as A~^. Therefore the 
non-Gaussian contributions for parallelogram configurations remain of the same order of 
magnitude as the Gaussian ones, with however, and again, an exception in the squeezed 
limit 9^0. 

The most interesting situations come when Y^ ^ 0. For these quadrilaterals, the 
Gaussian contribution vanishes and solely a non-Gaussian statistics can source the 
trispectrum. At large wavenumber (small angles), one has Y^ ^ 1 such that the 
normalised lower incomplete gamma function in Eq. f l95|) is close to unity: 

V^ L^ / - \ -l/(2x+2) 

Too(fci, fc2, fcs, ^4) ^ e^-^ [^c,e) fix)9iki, k,, k,, k,). (104) 

Under the scaling transformation ka — t- Afc^, the geometric factor, and thus Eq. (I104p . 
scales as 

g{\ki, Xk2, Afcg, Afc4) = X'''g{ki, ^2, ^3, ^4), (105) 

with 

p = 6 + ^. (106) 
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Figure 20. Trispectrum geometrical factor for tire kite quadrilaterals (represented on 
the left panel) as a function of the opening angle 9, and plotted for various values of a. 
The trispectrum is enhanced in the squeezed limit — > 0. The bump for Op = n — 2a 
corresponds to the parallelogram limit for which the unconnected part is no longer 
vanishing. 



As claimed, for NG strings, p is directly given by the power law of the the loop 
distribution |144j . Since this exponent is different from the one associated with 
parallelogram configurations it may actually be used to distinguish a trispectrum sourced 
by cosmic strings with the one generated by other non-Gaussian effects. 

In Fig. [201 "we have plotted the geometrical factor (7(^1,^2,^2,^4) for the kite 
quadrilaterals (represented in the same figure), as a function of their opening angles 9 
and a. As for the bispectrum, the trispectrum is enhanced on the squeezed configuration 
obtained when becomes small. In this limit Eq. ( 1971) can be expanded as 



8cos^( 



a] 



9 



kpe"' 



2 cos 2a.) 



X [2(1 + x)tan2(a) - 1 + 4^(1 -tan^ a 



-l/(2x+2) 



(107) 



and one recovers the mode dependency in k~^ while the amplitude is amplified as 9^^^^". 
As discussed in the previous section, the singular limit — )■ is never reached with a 
finite resolution beam. 



4-9. Comparison with data 

The cosmic string bispectrum and trispectrum associated with the flat polygonal 
configurations are the best suited to look for string signatures. However, it is not 
easy to compare with existing constraints as much of the literature focuses on particular 
models of primordial non-Gaussianity. For instance, in the local type of primordial non- 
Gaussianities, the parameter /nl characterises the primordial bispectrum and maximal 
amplitude occurs for squeezed triangle configurations, as it is the case for the cosmic 
strings j214] . However, as a result of the CMB transfer functions, a given value of /nl 
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corresponds to oscillating damped patterns of the CMB temperature bispectrum, which 
are completely different of the power laws we have found for the string bispectrum 
at small scales. The current bounds on /nl being precisely obtained from template 
matching procedures, they cannot be apphed to the strings [l2l I215[ 1216] . For this 
reason, the parameters used to quantify primordial non-Gaussianities are not well suited 
here, precisely because we expect the string non-Gaussianities to be non-primordial. An 
efficient approach would be to use a template matching procedure with the formulae 
derived in the previous sections. Another approach might be to estimate what values 
the primordial parameters, such as /nl and tnl, would assume if the non-Gaussianities 
were actually due to strings. Notice that asking such a question would be close to find 
the best amplitude of a sine function to fit a power law. However, since primordial 
non-Gaussianities are and will be tested in CMB data anyway, one could answer this 
question by performing a Fisher matrix analysis along the lines of Refs. [53l [55] . 

5. Conclusion and perspectives 

The results presented in this article were essentially concerned with Nambu-Goto type of 
cosmic strings, which is the simplest realisation of a one-dimensional spatially extended 
object. As a result, they should not be blindly extrapolated to other types of string, 
although, as argued in Sec. [21 some of them are expected to be generic. In particular, 
due to the scaling of the long strings, cosmic string loops do not infiuence significantly 
the CMB observables. Changing the intercommuting probability is expected to rescale 
some of the presented results )2171I218] . but in a way which remains to be quantified. 

In Sec. [HI we have briefiy reviewed the current understanding of the cosmological 
evolution of a string network by means of FLRW numerical simulations, which is a 
non-trivial problem even for NG strings. Observable predictions crucially depend on 
this step. Numerical simulations can be avoided by making some assumptions on the 
string distribution but at the expense of introducing unnecessary extra parameters. 
When approximate analytical models are then used to derive observable predictions, 
one should keep in mind that the results are as uncertain as the values assumed for 
the additional parameters. Provided one is interested in length scales not affected by 
gravitational back-reaction effects, all of the statistical properties of a NG string network 
in scaling depends only on one unknown physical parameter: the string energy density 
per unit length U, not mor affl - 

In this context. Sec. [4] discusses the non-Gaussian effects induced by a cosmological 
string network in the CMB temperature anisotropics. We have shown that string 
induced CMB fiuctuations have a negative skewness and a non-vanishing kurtosis. On 
a CMB temperature anisotropics map, these non-Gaussianities imprint characteristic 
signatures in a multifractal analysis as well as in the gradient magnitude, both being 
more significant at small angles. This property is recovered in Fourier space: the CMB 
angular power spectrum decays at most as 1/i, for the large multipoles £, and strings 

ft The expansion rate is supposed to be known. 
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become the dominant sources of primary fluctuations. The skewness appears to be 
the direct consequence of the breaking of the time reversal symmetry in an expanding 
universe, and implies the existence of a non-vanishing bispectrum. Using analytical 
approximations, tested and confirmed by numerical simulations, we then derived the 
expected bispectrum and trispectrum of string induced CMB temperature anisotropies 
for the large multipoles. Although the bispectrum decays not faster than £~^, the 
trispectrum multipole dependency is in (i~P, where p = 6 + l/(x + 1) and x is a small 
number related to the tangent vector correlator and the NG loop distribution. Due to the 
line-like CMB patterns induced by the strings, both the bispectrum and trispectrum are 
enhanced on all elongated triangle and quadrilateral configurations of the wavevectors. 
These ones may constitute the best configurations to look for a non-Gaussian string 
signal while being experimentally limited by finite beam resolution. Let us note that our 
expressions have been derived in the fiat sky approximation. String non-Gaussianities 
at small multipoles is still an open problem which could be dealt with full sky string 
maps |198l 12191 1220] . However, if, as the current constraints suggest, cosmic strings 
marginally contribute to the large scale CMB anisotropies, then they should show up 
at large multipoles in all of the above mentioned observables. This is precisely where 
the experimental efforts are directed. At very small angular scales, the difficulties will 
certainly be to separate the string signals from the astrophysical sources. Interestingly, 
the very soon accessible intermediate angles probed by the Planck satellite, and the 
other ground based telescopes, may not suffer from this problem and could be an open 
window on cosmic strings. 
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